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The azimuthal anisotropies of the collective transverse flow of hadrons are investigated in a large 
range of heavy-ion collision energy within the Parton-Hadron-String Dynamics (PHSD) microscopic 
transport approach which incorporates explicit partonic degrees of freedom in terms of strongly 
interacting quasiparticles (quarks and gluons) in line with an equation-of-state from lattice QCD 
as well as dynamical hadronization and hadronic dynamics in the final reaction phase. The experi- 
mentally observed increase of the elliptic flow V2 with bombarding energy is successfully described 
in terms of the PHSD approach in contrast to a variety of other kinetic models based on hadronic 
interactions. The analysis of higher-order harmonics «3 and U4 shows a similar tendency of grow- 
ing deviations between partonic and purely hadronic models with increasing bombarding energy. 
This signals that the excitation functions of azimuthal anisotropies provide a sensitive probe for the 
underling degrees of freedom excited in heavy-ion collisions. 
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Introduction. A few decades of experimental studies 
at the Schwcrioncn-Synchrotron (SIS), the Alternating 
Gradient Synchroton (AGS) and the Super Proton Syn- 
chroton (SPS) have shown that the physics of nuclear 
collisions at moderate relativistic energies is dominated 
by the nonequilibrium dynamics of hadronic resonance 
matter, i.e. the confined phase of QCD. The body of 
data extends and builds up the knowledge gained about 
dense hadronic matter, in particular, at the SPS/CERN. 
The SPS heavy-ion data have shown several signatures 
that hinted at the onset of a quark-gluon plasma (QGP) 
formation [1, 2]. With the Relativistic Heavy Ion Collider 
(RHIC) the center-of-mass energy could be increased by 
a factor of 10 relative to the SPS, and the experiments 
at RHIC assured that a new form of matter - well above 
the deconfinement transition point - was created in the 
laboratory. 

Indeed, the discovery of a large azimuthal anisotropic 
flow of hadrons at RHIC provides a conclusive evi- 
dence for the creation of dense partonic matter in ultra- 
relativistic nucleus-nucleus collisions. The strongly inter- 
acting medium in the collision zone can be expected to 
achieve a local equilibrium and exhibit an approximately 
hydrodynamic flow [3-5]. The momentum anisotropy is 
generated due to pressure gradients in a collective expan- 
sion of an initial geometry of an "almond-shaped" colli- 
sion zone produced in noncentral collisions [3, 4]. The 
pressure gradients translate early stage coordinate space 
asymmetry to final-state momentum space anisotropy [6] . 
The picture thus emerges that the medium created in 
ultra-relativistic collisions for a couple of fm/c interacts 
more strongly than hadron resonance matter and exhibits 
collective properties that resemble those of a liquid of a 
very low shear viscosity rj to the entropy density s ratio, 
r]/s, close to a nearly perfect fluid [7-9]. 



An experimental manifestation of this collective flow is 
the anisotropic emission of particles in the plane trans- 
verse to the beam direction. A quark number scaling 
of the elliptic flow proposed in Ref. [10] was observed 
at RHIC for a broad range of particle species, collision 
centralities, and transverse kinetic energy, presumably 
interpreted as due to the development of substantial col- 
lectivity in the early partonic phase [11]. 

It was shown that higher-order anisotropy harmonics, 
in particular the hexadecupole moment V4, can provide 
a more sensitive constraint on the magnitude of 77/s and 
the freeze-out dynamics, and the ratio W4/(w2)^ might 
indicate whether a full local equilibrium is achieved in 
the QGP [12]. Recently, the importance of the trian- 
gular flow W3, which originates from fluctuations in the 
initial collision geometry, has been pointed out [13, 14]. 
The participant triangularity characterizes the triangu- 
lar anisotropy of the initial nuclear overlap geometry and 
arises from event-by-event fluctuations in the participant- 
nucleon collision points and corresponds to a large third 
Fourier component in two-particle azimuthal correlations 
at large pseudo-rapidity separation A?/. This fact sug- 
gests a significant contribution of the triangular flow to 
the ridge phenomenon and broad away-side structures 
observed in the RHIC data [13]. 

A large number of anisotropic flow measurements 
have been performed by many experimental groups at 
SIS, AGS, SPS and RHIC energies over the past 20 
years. Very recently the azimuthal asymmetry has also 
been measured at the Large Hadron Collider (LHC) at 
CERN [15]. The Beam Energy Scan (BES) program 
proposed at RHIC [16] covers the energy interval from 
y/sNN = 200 GeV, where partonic degrees of freedom 
(DOF) play a decisive role, down to the AGS energy 
y^SNN ~ 5 GeV, where most experimental data can be 
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described successfully in terms of hadronic DOF Lower- 
ing the collision energy and studying the energy depen- 
dence of an anisotropic flow allows one to search for the 
onset of the transition to a phase with partonic DOF at 
an early stage of the collision as well as possibly iden- 
tify the location of the expected critical end-point that 
terminates the first order phase transition at high quark- 
chemical potential [11, 17]. 

This work aims to study excitation functions for dif- 
ferent harmonics of the charged particle anisotropy in 
momentum space in a wide collision energy range, i.e. 
from the AGS to the top RHIC energy regime. We want 
to clarify how the interplay of quark and hadron DOF 
is changed with increasing bombarding energy. In this 
study we investigate the excitation function of different 
flow coefficients. Our analysis of the STAR/PHENIX 
RHIC data ~ based on recent results of the BES pro- 
gram - will be performed within the PHSD transport ap- 
proach [18] that includes explicit partonic DOF as well as 
a dynamic hadronization scheme for the transition from 
partonic to hadronic DOF and vice versa. 

The PHSD approach. The dynamics of partons, 
hadrons and strings in relativistic nucleus-nucleus col- 
lisions is analyzed here within the Parton-Hadron-String 
Dynamics approach [18] . In this transport approach the 
partonic dynamics is based on Kadanoff-Baym equations 
for Green functions with self-energies from the Dynami- 
cal QuasiParticle Model (DQPM) [19, 20] which describes 
QCD properties in terms of "resummed" single-particle 
Green functions. In Ref. [21], the actual three DQPM 
parameters for the temperature-dependent effective cou- 
pling were fitted to the recent lattice QCD results of 
Ref. [22]. The latter lead to a critical temperature Tc « 
160 MeV which corresponds to a critical energy density 
of ec ~ 0.5 GeV/fm'^. In PHSD the parton spectral func- 
tions Pj (j = 5, q, g) are no longer ^-functions in the in- 
variant mass squared as in conventional cascade or trans- 
port models but depend on the parton mass and width 
parameters which were fixed by fitting the lattice QCD 
resuhs from Ref. [22]. We recall that the DQPM al- 
lows one to extract a potential energy density Vp from 
the space-like part of the energy-momentum tensor as 
a function of the scalar parton density p^. Derivatives 
of Vp with respect to ps then define a scalar mean-field 
potential Us{ps) which enters into the equation of mo- 
tion for the dynamic partonic quasiparticles. Further- 
more, a two-body interaction strength can be extracted 
from the DQPM as well from the quasiparticle width 
in line with Ref. [9]. The transition from partonic to 
hadronic DOF (and vice versa) is described by covariant 
transition rates for the fusion of quark- antiquark pairs 
or three quarks (antiquarks) , respectively, obeying flavor 
current-conservation, color neutrality as well as energy- 
momentum conservation [18, 21]. Since the dynami- 
cal quarks and antiquarks become very massive close to 
the phase transition, the formed resonant "prehadronic" 
color-dipole states {qq or qqq) are of high invariant mass, 
too, and sequentially decay to the ground-state meson 



and baryon octets increasing the total entropy. 

On the hadronic side PHSD includes explicitly the 
baryon octet and decouplet, the 0~- and 1 "-meson 
nonets as well as selected higher resonances as in 
the Hadron-String-Dynamics (HSD) approach [23, 24]. 
Hadrons of higher masses (> 1.5 GeV in case of baryons 
and > 1.3 GeV for mesons) are treated as "strings" (color 
dipolcs) that decay to the known (low-mass) hadrons, ac- 
cording to the JETSET algorithm [25]. Note that PHSD 
and HSD merge at low energy density, in particular below 
the critical energy density Cc « 0.5 GeV/fm'^. 

The PHSD approach was applied to nucleus-nucleus 
collisions from -y/sTvw ^ 5 to 200 GeV in Refs. [18, 21] 
in order to explore the space-time regions of "partonic 
matter" . It was found that even central collisions at the 
top-SPS energy of ^/sNN = 17.3 GeV show a large frac- 
tion of nonpartonic, i.e., hadronic or stringlike matter, 
which can be viewed as a hadronic corona. This flnding 
implies that neither hadronic nor only partonic "mod- 
els" can be employed to extract physical conclusions in 
comparing model results with data. All these previous 
flndings provide promising perspectives to use PHSD in 
the whole range from about ^/snTj = 5 to 200 GeV for 
a systematic study of azimuthal asymmetries of hadrons 
produced in relativistic nucleus- nucleus collisions. 

Calculational results and comparison to data. The 
anisotropy in the azimuthal angle ip is usually char- 
acterized by the even order Fourier coefficients w„ = 
{exp{ I n{4> — ^'flp))), = 2, 4, since for a smooth an- 
gular proflle the odd harmonics become equal to zero. 
As noted above, \E'flp is the azimuth of the reaction 
plane and the brackets denote averaging over particles 
and events. In particular, for the widely used second- 
order coefficient, denoted as an elliptic flow, we have 

V2 = {cos{2^ - 2^np)) = ( ) ' (1) 

\Px + Py I 

where px and Py are the x and y components of the par- 
ticle momenta. This coefficient can be considered as a 
function of ccntrality, pseudorapidity 77, and/or trans- 
verse momentum pT. We note that the reaction plane in 
PHSD is given by the {x — z) plane with the z axis in 
the beam direction. Integrated W3 and Vi coefficients are 
calculated by the two-particle correlation method in line 
with Ref. [26]: 

{cosijiipi - nip2))) = (vl) + Sn , (2) 

which is based on the assumption that nonfiow contribu- 
tions 6n are small and account for event-by-event fluctu- 
ations of the event plane. In Fig. 1 the experimental V2 
excitation function in the transient energy range is com- 
pared to the results from the PHSD calculations; HSD 
model results are given as well for reference. We note 
that the centrality selection and acceptance are the same 
for the data and models. 

We recall that the HSD model has been very success- 
ful in describing heavy-ion spectra and rapidity distribu- 
tions from SIS to SPS energies. A detailed comparison 
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FIG. 1: Average elliptic flow V2 of charged particles at 
midrapidity for two centrality selections calculated within the 
PHSD (solid curves) and HSD (dashed curves) approaches. 
The V2 STAR data compilation for minimal bias collisions 
are taken from Ref. [27] (stars) and the preliminary PHENIX 
data [28] are plotted by solid circles. 



of HSD results with respect to a large experimental data 
set was reported in Ref. [29] for central Au+Au (Pb+Pb) 
collisions from AGS to top SPS energies. Indeed, as 
shown in Fig. 1 (dashed lines), HSD is in good agree- 
ment with experiment for both data sets at the lower 
edge (y^Sjviv ~ 10 GeV) but predicts an approximately 
energy-independent flow V2 at larger energies and, there- 
fore, does not match the experimental observations. This 
behavior is in quite close agreement with another inde- 
pendent hadronic model, the UrQMD (Ultra relativistic 
Quantum Molecular Dynamics) [30] (cf. with [27]). 

From the above comparison one may conclude that the 
rise of V2 with bombarding energy is not due to hadronic 
interactions and models with partonic DOF have to be 
addressed. Indeed, the PHSD approach incorporates the 
parton medium effects in line with a IQCD equation of 
state, as discussed above, and also includes a dynamic 
hadronization scheme based on covariant transition rates. 
It is seen from Fig. 1 that PHSD performs better: The 
elliptic flow V2 from PHSD (solid curve) is fairly in line 
with the data from the STAR and PHENIX collabora- 
tions and clearly shows the growth of V2 with the bom- 
barding energy. 

Since partonic DOF come into play, it is interesting 
to compare with another parton-hadron model, i.e. the 
AMPT (A Muhi Phase Transport) model [31]. This 
model is based on a perturbative QCD description of par- 
tonic interactions, including the production of multiple 
minijet partons according to the number of binary initial 
collisions. As shown in Ref. [27], the AMPT model pre- 
dicts an approximately constant V2 with i/snn similar 
to the hadronic models HSD and UrQMD; however, the 
V2 values match the experimental data at the top RHIC 
energy. This discrepancy is due to a pQCD description of 
the partonic phase in AMPT where the minijet partons 
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FIG. 2: Average anisotropic flows vs and V4 of charged par- 
ticles at mid-pseudorapidity for minimum bias collisions of 
Au + Au calculated within the PHSD (solid lines) and HSD 
(dashed lines) models. 



are treated as massless and their potentials are disre- 
garded when they undergo scattering. Note that PHSD 
and AMPT (with the additional strong hadron melting 
assumption in AMPT) practically give the same ellip- 
tic flow at the top RHIC energy of y^SNN = 200 GeV. 
We note that the PHSD model includes more realistic 
properties of dynamical quasiparticles especially in the 
vicinity of the critical energy density. Furthermore, the 
quark-gluon transport in PHSD naturally passes on to 
the (hadronic) HSD model at lower y^SNN- 

In Fig. 2 we display the PHSD and HSD results for 
the anisotropic flows V3 and U4 of charged particles at 
midpseudorapidity for Au + Au collisions from ^/snn = 
5 to 200 GeV. The triangular flow increases with y^s^N 
having negative values for ^/snn ~ 10 GeV. The pure 
hadronic model HSD gives V3 ~ for ^/snn ~ (20- 
30) GeV. Accordingly, the results from PHSD (solid red 
lines) are systematically larger than those from HSD 
(dashed blue lines). Unfortunately, our statistics is not 
high enough to allow for more precise conclusions. The 
hexadecupole flow W4 stays almost constant in the con- 
sidered energy range; here PHSD gives slightly higher 
values than HSD. 

The V2 increase is clarified in Fig. 3 where the partonic 
fraction of the energy density at mid-pseudorapidity with 
respect to the total energy density in the same pseu- 
dorapidity interval is shown. We recall that the repul- 
sive scalar mean field potential Us{ps) for partons in the 
PHSD model leads to an increase of the flow V2 as com- 
pared to that for HSD or PHSD calculations without par- 
tonic mean fields [21]. As follows from Fig. 3, the energy 
fraction of the partons substantially grows with increas- 
ing bombarding energy while the duration of the partonic 
phase is roughly the same. Accordingly, the increasing in- 
fluence of the repulsive partonic mean-field Us{ps) leads 
to an increase of the flow V2 with bombarding energy. 
We point out that the increase of 1)2 in PHSD relative to 
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FIG. 3: The evolution of the parton fraction of the total 
energy density at the mid-pseudorapidity for different collision 
energies. 

HSD is also partly due to the higher interaction rates in 
the partonic medium because of a lower ratio oi rj/s for 
partonic degrees of freedom at energy densities above the 
critical energy density than for hadronic media below the 
critical energy density [32, 33]. The relative increase in 
and V4 in PHSD essentially is due to the higher partonic 
interaction rate and, thus, to a lower ratio rj/s in the 
partonic medium, which is mandatory to convert initial 
spacial anisotropics to final anisotropics in momentum 
space [34]. 

Conclusions. The anisotropic flows - elliptic f2, tri- 



angular U3, and hexadecupole Vi - are reasonably de- 
scribed within the PHSD model in the whole transient 
energy range naturally connecting the hadronic processes 
at moderate bombarding energies with ultrarelativistic 
collisions at RHIC energies where the quark-gluon DOF 
become dominant due to a growing number of partons. 
The smooth growth of the elliptic flow with collision 
energy demonstrates the increasing importance of par- 
tonic DOF This feature is reproduced by neither explicit 
hadronic kinetic models like HSD or UrQMD nor the 
AMPT model treating the partonic phase on the ba- 
sis of pQCD with massless partons and a noninteracting 
equation-of-state for the partons. Further signatures of 
the transverse collective flow, the higher-order harmonics 
of the transverse anisotropy and U4 change only weakly 
from y/SNN w 7 GeV to the top RHIC energy ^/SNN = 
200 GeV, roughly in agreement with preliminary exper- 
imental data. Certainly, new measurements within the 
BES program at RHIC, especially for higher-order har- 
monics, will further constrain the partonic dynamics. 
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